Effects of Bedrock Landsliding on Cosmogenically Determined Erosion Rates 
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The successful quantification of long-term erosion rates underpins our understanding of landscape formation, the topographic 
evolution of mountain ranges, and the mass balance within active orogens. The measurement of in situ-produced cosmogenic 
radionuclides (CRNs) in fluvial and alluvial sediments is perhaps the method with the greatest ability to provide such long- 
term erosion rates. In active orogens, however, deep-seated bedrock landsliding is an important erosional process, the effect 
of which on CRN-derived erosion rates is largely unquantified. We present a numerical simulation of cosmogenic nuclide 
production and distribution in landslide-dominated catchments to address the effect of bedrock landsliding on cosmogenic 
erosion rates in actively eroding landscapes. Results of the simulation indicate that the temporal stability of erosion rates 
determined from CRN concentrations in sediment decreases with increased ratios of landsliding to sediment detachment rates 
within a given catchment area, and that larger catchment areas must be sampled with increased frequency of landsliding in order 
to accurately evaluate long-term erosion rates. In addition, results of this simulation suggest that sediment sampling for CRNs is 
the appropriate method for determining long-term erosion rates in regions dominated by mass-wasting processes, while bedrock 
surface sampling for CRNs is generally an ineffective means of determining long-term erosion rates. Response times of CRN 
concentrations to changes in erosion rate indicate that climatically driven cycles of erosion may be detected relatively quickly 
after such changes occur, but that complete equilibration of CRN concentrations to new erosional conditions may take tens of 
thousands of years Simulation results of CRN erosion rates are compared with a new, rich dataset of CRN concentrations from 
the Nepalese Himalaya, supporting conclusions drawn from the simulation. 


1. Introduction 

The successful quantification of long-term ero- 
sion rates underpins our understanding of landscape 
formation, the topographic evolution of mountain 
ranges, and the mass balance within active orogens 
[1]. Observed changes in long-term erosion rates are 
often considered proxies for changes in the climatic or 
tectonic boundary conditions that control landscape 
evolution. As the importance of long-term erosion 
rates to unraveling significant geomorphic and tec- 
tonic problems has become clear, a widely applica- 
ble means of measuring long-term surface denudation 
rates has been sought. The measurement of in situ- 
produced cosmogenic radionuclides (CRNs) in fluvial 
and alluvial sediments has been shown to yield spa- 
tially averaged erosion rates, and has become perhaps 
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the method with the greatest applicability in measur- 
ing erosion rates over 1 0 3 — 1 0 5 years and across a 
wide variety of landscapes and erosional processes 
|2— 6]. Most existing studies that utilize CRN -derived 
erosion rates have focused on regions with fairly 
spatially and temporally homogeneous erosion rates. 
Here we present a simulation to explore the effects of 
bedrock landsliding on cosmogenic erosion rates, and 
potential for exploiting cosmogenic nuclides to mea- 
sure erosion rates in rapidly eroding, active orogens. 

Large bedrock landslides incise to depths greater 
than one or more attenuation lengths of cosmic rays, 
thus mobilizing sediments with little or no cosmo- 
genic nuclide abundance (Fig. 1). For example, based 
on empirical relationships of landslide depth to area 
17], a landslide with a radius of just 10 meters will in- 
cise to ~ 100 cm, below one attenuation depth for 
spallogenic nuclide production. An extraordinarily 
large slide may incise to a maximum depth of sev- 
eral tens of meters, below one attenuation depth for 
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Figure 1. Cartoon showing the effects of landslides on surface cosmogenic nuclide concentrations. A) Homogeneous CRN surface concen- 
tration and sediment volume and concentration during steady-state erosion. B) Heterogeneous surface CRN concentration and heterogeneous 
sediment volume and CRN concentration during landslide-dominated erosion. 


muogenic nuclide production. In catchments where 
deep landslides have recently occurred, the addition 
of nuclide-poor landslide detritus to the fluvial system 
will dilute the CRN concentration in the fluvial sed- 
iment, yielding apparently higher erosion rates. For 
example, samples from low-order fluvial catchments 
in the Nepalese Himalaya yield ‘erosion rates’ de- 
termined from cosmogenic nuclides that range from 
0.01 to 0.5 mm/yr, even along the same ridge crest 
(Heimsath, in prep.). Similar results are reported from 
the San Bernardino Mountains in southern California, 
where ‘erosion rates’ along a landslide-dominated es- 
carpment vary from < 0.3 to > 2.7 mm/yr, while rates 
from an adjacent region dominated by sediment de- 
tachment are on the order of hundredths of mm/yr 
18]. Whereas such variability may be expected for 
data sets that are focused on small, steep catchments 
in actively deforming mountain belts, it does not ex- 
plain how in situ-produced CRNs may be utilized to 
understand the rates of erosion in regions affected by 
both sediment detachment and landsliding processes. 

Here we address the effects of bedrock landslid- 
ing on CRN-derived erosion rates from bedrock and 


fluvial samples in an attempt to answer several basic 
questions: As the rate of landsliding increases, how 
are CRN concentrations in fluvial sediment affected? 
Can reliable CRN erosion rates be derived from flu- 
vial sediment or bedrock samples when landsliding 
is the dominant erosional process? Over what spatial 
scales do fluvial systems integrate the effects of land- 
sliding? What timescales are required for CRN con- 
centrations to respond to changes in erosion rates? We 
present the results of a numerical simulation of cos- 
mogenic nuclide production and erosional removal in 
landslide-dominated catchments to assess these ques- 
tions. A series of simulations with varying sediment 
detachment and landslide erosion rates are used to 
create statistical populations of CRN-derived erosion 
rates for both ‘sediment’ and ‘bedrock’ samples for a 
theoretical landscape. Finally, simulated distributions 
of CRN-derived ages are compared to a new', rich data 
set of CRN-derived erosion rates from the Nepalese 
Himalaya. 
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2. Numerical Simulation 

The numerical simulation is based on actual digital 
elevation data, and simulates the production of cos- 
mogenic nuclides at each model cell, the removal of 
material through sediment detachment, the removal 
of material by landsliding, and the radioactive de- 
cay of cosmogenic nuclides. A Geographic Informa- 
tion System (GIS; in this case Arclnfo) is used for 
the backbone of the simluation. Model initialization, 
data assimilation, and data output are all controlled 
through the GIS using Arc Macro Language (AML). 
Computationally intensive portions of the model are 
passed from the GIS to customized Perl modules for 
computational efficiency. The three main functional- 
ities of the simulation are described in greater detail 
below. 


2.1. Cosmogenic Nuclide Production 

Prior to running the landslide simulation model, 
cosmogenic nuclide production rates must be calcu- 
lated for each cell in the model. Calculation of these 
rates begins with a geo-referenced digital elevation 
model (DEM) of the study area of interest. For each 
cell in the digital elevation model, cosmogenic pro- 
duction scaling factors are calculated within the GIS 
based on cell altitude and latitude following [9,10]. 
Further corrections to cosmogenic production are ap- 
plied by calculating the topographic shielding at each 
point in the DEM. For each cell in the DEM, the ver- 
tical angle to to every other cell is calculated. These 
values are binned into 5° radial bins, and the maxi- 
mum vertical angle in each bin is used to approximate 
the horizon angle for that bin. The topographic shield- 
ing factor for each bin is derived from the horizon 
angle using a published methodology [11]. The al- 
titude, latitude, and topographic shielding factors are 
combined within the GIS system to produce an output 
array of cosmogenic production scaling factors. This 
array is then multiplied by the high-latitude, sea-level 
production rate of the cosmogenic nuclide of interest 
to create an array of cosmogenic production rates. In 
this case, we have chosen to model J0 Be and selected 
a production rate of 5.3 atoms/g/yr [12]. This array of 
10 Be production rates is preserved for use through the 
rest of the model run. 


22. Model Initialization 

Model initialization consists of two separate ac- 
tions: i) preparing the model for data gathering and 
assimilation and ii) calculating an initial surface cos- 
mogenic nuclide concentration to start the landslide 
model. The first of these two tasks is the most time- 
intensive and must be performed separately for each 
DEM on which this model is run. Using the hydro- 
logic functions available in a GIS, watersheds within 
the model area are delineated. A variety of first- 
through highest-order watersheds are selected and 
saved for later use in data analysis. Additionally, 100 
random points are generated across the model space 
at which to track cosmogenic nuclide concentrations 
in bedrock and erosional removal of material. After 
the data gathering initialization steps are complete, an 
input cosmogenic surface nuclide concentration array 
must be calculated for input to the landslide portion of 
the model. We have chosen to input a surface concen- 
tration grid that represents the steady-state concentra- 
tion of 10 Be at the sediment detachment rate specified 
for the model run. The sediment detachment rate is 
limited by the rate at which rock can be converted to 
soil or regolith, and is taken to be less than 0.3 mm/yr 
[13]. The surface concentration ( N ) is calculated for 
each point in the mode! following this equation: 


N 10 Be = 


P 10 Be 

(Aio Be -f E/A) 


(1 _ e ~(^Be + E / A '> t ) (]) 


Where jV 10 Be and P 10 Be are the concentration 
(atoms/g) and production rate (atoms/g/yr) of 1() Be, 
respectively; Aio Bc is the decay constant of 10 Be 
(yr _1 ); E is the erosion rate (here equivalent to the 
sediment detachment rate), in g/cm 2 /yr; A is the neu- 
tron attenuation length in rock, in g/cm 2 ; and t is time 
(yr). This initial concentration array is saved for input 
into future model runs. 


2 _3. Landslide Simulation 

The landslide simulation portion of the model takes 
as input the surface cosmogenic nuclide concentra- 
tion array calculated in the previous step, and the cos- 
mogenic production rate array calculated in the first 
model step. This portion of the model is iterated. 
At the end of each iteration, two arrays are output, a 
depth array which contains the sum of all sediment 
removed by erosion, due to both sediment detach- 
ment and landsliding, and a surface concentration ar- 
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ray that contains the surface cosmogenic nuclide con- 
centration at each model cell after erosional removal 
of material, cosmogenic ingrowth and radioactive de- 
cay. The individual steps are detailed below. 

23.1. Cosmogenic ingrowth and decay 

The first step in the iteration of the landslide model 
is to add cosmogenic nuclide ingrowth for one model 
time step (t) to the input surface cosmogenic nuclide 
concentration array. Given, for each cell in the model, 
an initial surface cosmogenic nuclide concentration, 
./V t 10 Be, and a surface cosmogenic nuclide production 
rate, P 10 Be, a resultant surface cosmogenic nuclide 
concentration, AfJ°Be, produced by cosmic ray bom- 
bardment and removal by radioactive decay can be 
calculated by: 

AT r 10 Be = (iV/ 0 Be + P 10 Be • t)e~ Xl0 B «* (2) 

The new surface concentration, 7V r 10 Be, replaces 
the initial value in the surface concentration array. 

2.3.2. Sediment detachment 

For each model run, a sediment detachment rate, 
E s is specified. This rate represents the spatially ho- 
mogeneous erosional removal of material from the 
land surface, limited by weathering and soil produc- 
tion processes (e.g. grussification in granitic ter- 
ranes). For each model time step (t), a depth equal 
to E s x t is removed from the landscape and added to 
the depth grid. This material is always removed from 
the upper surface of the topography, and therefore has 
the highest concentration of CRNs. 

2.3.3. Landslides 

After removal of material by sediment detach- 
ment, the model is populated with landslides. Land- 
slides are assumed to obey a power-law frequency- 
magnitude relationship [7,14,151. Based on this as- 
sumption, populations of landslides in the model can 
be derived from four parameters: /3, the power-law 
exponent for landslide frequency-magnitude relation- 
ship, A m i n , the minimum landslide area considered 
in the model, A max , the maximum landslide area con- 
sidered in the model, and Ei s , the average rate of ero- 
sion by landsliding over the model area. Although 
short-term landslide erosion rates fluctuate due to the 
episodicity of landslides, the average rate of ero- 
sion by landsliding is produced through the power- 
law frequency-magnitude relationship over time (Fig. 



Figure 2. Comparison of output model landslide erosion rates with 
prescribed landslide erosion rates, averaged over 10,000 year time 
steps. Mean erosion rates for the 1 00,000 year model run are shown 
at the right of the figure. Short term variations in landslide erosion 
rate are due to random distribution of slide size; long term landslide 
erosion rates return input erosion rates ±3%. 


2). Complete derivations of landslide frequency- 
magnitude distributions in the model space are given 
in Appendix A. For each model timestep, a land- 
slide distribution is generated, and the landslides are 
randomly distributed over the model space. The to- 
tal amount of material removed by landsliding from 
each model cell during the timestep is then calculated. 
The total depth of material removed by landsliding is 
added to the depth of material removed by sediment 
detachment (see Appendix A). 

2.3.4. Surface concentration 

Once the depth array is tabulated, the surface cos- 
mogenic nuclide concentration is recalculated to re- 
flect the depth of material removed from each model 
cell. Assuming, for simplification, a constant rock 
density ( p ) in the model of 2.65 g/cm 8 , and a cos- 
mogenic ray attenuation length in this rock (A) of 150 
cm 2 /g, a final surface cosmogenic nuclide concentra- 
tion at the end of the timestep, Nj° Be. can be cal- 
culated for each cell based on the cosmogenic con- 
centration following ingrowth, A r r 10 Be. above, and the 
depth of material. D. removed from the cell during the 
timestep: 

A r y°Be = A r r 10 Be ■ e~ Dip/A '> (3) 

The value Nj°Be is stored for each cell in the surface 
concentration array. At this point, this array is saved 
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for data extraction, as described in the next section, 
and then re-opened as the initial surface concentration 
grid for the next iteration step of the landslide portion 
of the model. 


2.4. Data Extraction 

At a user-specified sampling interval, the depth ar- 
ray and surface concentration array are passed back 
to the GIS, and data are extracted from each array to 
simulate two potential cosmogenic nuclide sampling 
methods: surface exposure age dating and stream sed- 
iment sampling. 

2.4.1. Surface exposure age dating 

To simulate the effects of landsliding on surface 
age exposure dating, at each sampling interval, the 
depth of material removed from, and the final sur- 
face concentration at, each of 100 randomly gener- 
ated points across the model are recorded. From these 
data, a volumetric erosion rate ( E v ) for this sampling 
interval is calculated by dividing the depth of material 
removed ( D ) by the model time step (f; recall that D 
includes both continual removal of rock by sediment 
detachment and episodic landsliding). Additionally, a 
cosmogenic erosion rate ( E c ) is calculated, for each 
point using the standard assumption of steady-state 
erosion [16]: 


E c 


A ( P 10 Be 
p \N { W Be Al ° 


(4) 


Both volumetric and cosmogenic surface-exposure- 
age erosion rates, along with the depth and surface 
concentration information, are stored in a database 
file. Data can be extracted after the completion of 
the model run to analyze the variation in volumetric 
and cosmogenic erosion rates at a given point through 
time, or data from all points can be assimilated and a 
probability density function (PDF) for either erosion 
rate can be calculated. 


2.4.2. Stream sediment sampling 

In addition to surface age exposure dating, cosmo- 
genic nuclide concentrations in stream sediments can 
be used to estimate average upstream erosion rates. 
To test the effects of landsliding on such estimates of 
erosion rates, the model also extracts data by water- 
sheds. At each sampling interval, the following pro- 
cedures are carried out for each watershed of inter- 
est. A volumetric erosion rate (E v ) is calculated by 


summing the total depth of material removed from the 
watershed and dividing it by the watershed area mul- 
tiplied by the model time step. Second, a cosmogenic 
erosion rate ( E c ) for the watershed is calculated. To 
calculate this erosion rate, first the concentration of 
the cosmogenic nuclide of interest in the eroded sedi- 
ment must be determined. Using the depth array and 
surface concentration array, the average nuclide con- 
centration in the eroded material at each model cell, 
A r l, 0 Be, can be calculated as follows: 

AT^Be = AT/ 10 Be • A (l - e - D(p/A) ) (5) 

The concentration of cosmogenic nuclides in the sed- 
iment removed from the watershed, AT| 0 Be, then, is: 

^ 0 Be = ^A“Be/^P (6) 

The erosion rate derived from the cosmogenic nuclide 
concentration in the material removed from the water- 
shed is then calculated as in Eqn. 4, where P 10 Be 
would represent the watershed-averaged production 
rate (derived from the GIS-based production calcula- 
tions), and A^°Be would be replaced with Ng°Be. 
Both volumetric and cosmogenic erosion rates are 
recorded to database files, and saved for later anal- 
ysis. 

2.43. Simplifications and assumptions 

It should be noted that the model outlined above 
makes several simplifying assumptions regarding sed- 
iment production and transport. With regards to sedi- 
ment detachment and transport, this is assumed to oc- 
cur at a homogeneous rate across the entire model, 
although in reality the rate at which these processes 
act on the scale of our model are controlled by lo- 
cal slope and lithology. Landslides are also randomly 
placed across the landscape, with no considerations 
for hillslope and aspect. Second, the model contains 
no provision for sediment storage. All material de- 
rived from landsliding is assumed to pass through 
the model within the timestep in which the landslide 
occurred (100 years for the models discussed here). 
Third, the model is intended to produce a population 
of cosmogenic and volumetric erosion rates for statis- 
tical analysis. Erosion and landslides occur through 
‘time’ to produce a variety of surface CRN concen- 
trations that could potentially be sampled; however, 
this model is not a landscape evolution model. Dur- 
ing the course of the model run, the model landscape 
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Table 1 

Statistics for basins tracked in the model. 


Basin Order 

Number of 
Basins 

Mean Area 
(km 2 ) 

First 

30 

0.1 

Second 

6 

0.6 

Third 

5 

2.3 

Fourth 

4 

8.4 

Fifth 

1 

69.0 


surface does not evolve, and shielding effects or ab- 
solute elevation changes that in reality would alter 
CRN production rates are not considered. Finally, the 
role of muogenic production is not explicitly consid- 
ered in this simulation. The greater attenuation length 
of muons would likely have the effect of moderating 
to some extent the surface nuclide concentration in 
the wake of small- to moderate-sized landslides, but 
large landslides, which carry the majority of the ero- 
sional load, will incise more deeply than the attenua- 
tion depth of muons, and thus would have a negligible 
effect on the model. 

3. Model Results 

To examine the model, we simulated the effects 
of landsliding on cosmogenic nuclide equilibrium, 
and associated CRN-derived and volumetrically cal- 
culated erosion rates for the San Antonio Creek catch- 
ment, located in the eastern San Gabriel Mountains of 
southern California (Fig. 3). It is a small (~ 70 km 2 ), 
mountainous catchment, selected in part because of 
the availability of high-quality digital elevation model 
(DEM) data over the region (30-meter resolution), 
and in part because a significant amount of work ex- 
ists describing the geomorphology and neotectonics 
of the region. Low -temperature thermochronologic 
data indicate that the eastern San Gabriel Mountains 
are being exhumed at a rate of ~ 0.3 — 1 mm/yr f 17- 
19], while geomorphic and geologic studies indicate 
that landsliding is a prevalent mechanism of erosion 
in this watershed [20,15]. In fact, both a landslide 
frequency-magnitude exponential scaling factor (/?) 
and a long-term average erosion rate have been de- 
termined for this region [15]. Interestingly, in the San 
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Figure 3. Shaded relief map of the San Antonio Creek watershed. 
San Gabriel Mountains. California. The full watershed is outlined 
with a heavy black line, while representative smaller order water- 
sheds are indicated with a lighter black line. Random ‘bedrock’ 
sampling points (discussed in the text) are shown as black dots. 


Gabriel Mountains, J3 = 1.1, similar to results from 
the Southern Alps of New Zealand [7]. When 0 is less 
than 1.5, large, but infrequent landslides dominate the 
overall sediment flux from a catchment. As such, this 
watershed potentially provides a natural laboratory to 
study the effects of landsliding on CRN-derived ero- 
sion rates, and for comparison with and calibration of 
the numerical model. 

Within the San Gabriel Mountains, San Antonio 
Creek is a fifth-order stream. The San Antonio Creek 
watershed was divided into 46 sub-basins that were 
tracked as part of this simulation, from which mod- 
eled sediment-derived CRN erosion rates were calcu- 
lated (examples shown in Fig. 3; basin statistics are 
listed in Table 1). In addition to the sub-basins, 100 
points were randomly distributed across the model 
space to serve as simulated bedrock CRN sampling 
localities (Fig. 3). For San Antonio Creek, sedi- 
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Figure 4. Example of a box -and- whisker plot to graphically display 
the statistical distribution of a non-Gaussian data set. A probability 
density function (PDF) and boxplot are shown for actual model 
results of cosmogenic erosion rates from first-order basins with a 
sediment detachment rate of 0.1 mm/yr and a total erosion rate of 
10 mm/yr. The distribution of erosion rates is plotted across the 
bottom of the graph as a series of vertical black hash marks. 


ment detachment rates of 0.01, 0.1 and 0.3 mm/yr 
were selected, spanning the range of observed sedi- 
ment detachment rates in the Transverse Ranges [8] 
to the most rapid known rates of rock-to-regolith con- 
version [13]. For each of the three sediment detach- 
ment rates, simulations of total erosion rates of 1, 
5, and 10 mm/yr were run. The rate of erosion by 
landsliding in each simulation run is the difference 
between the total erosion rate and the sediment de- 
tachment rate. Once the simulation was equilibrated 
to the imposed erosion rate, the simulation ran for 
100,000 yrs in 100 yr timesteps. Simulated cosmo- 
genic nuclide concentrations were sampled from the 
46 tracked catchments and the 100 bedrock locations 
every 1000 years, yielding 100 data points per analy- 
sis element (either catchment or bedrock location). 

Results from the simulation are summarized and 
displayed as a box-and-whisker plot. Such a plot al- 


lows a fairly straight-forward visual presentation of 
the statistical distribution of a data set, and is a valu- 
able graphical method to compare the distribution of 
two separate data sets (Fig. 4). 

3.1. Simulated sediment erosion rates 

Catchment-wide erosion rates from the 46 catch- 
ments tracked in the model allow comparison of 
the statistical distributions between CRN-derived and 
volumetrically averaged erosion rates for each of the 
nine simulations, and illustrate variations in the statis- 
tical distributions within each simulation as a function 
of catchment order (Fig. 5). 

Several conclusions can be drawn from these data. 
First, with increasing proportion of sediment detach- 
ment, the CRN-determined erosion rate at any catch- 
ment scale more closely reflects the volumetric ero- 
sion rate for any given total combined erosion rate 
(Fig. 6). This is not unexpected, as sediment de- 
tachment is modeled as a uniformly continuous pro- 
cess in the simulation, such that increased ratios of 
sediment detachment-to-landsliding will result in a 
greater contribution from a steady-state process to 
the overall erosional volume. (It is worth noting 
that under this formulation, for any given total ero- 
sion rate, an increase in sediment detachment rate re- 
sults in a decrease in the rate of erosion due to land- 
sliding. This, increasing the rate of sediment de- 
tachment CRN-derived erosion rates that are closer 
to the volumetric erosion rates, but further from the 
total erosion rate due the reduction in frequency of 
landslides; Fig. 6). Additionally, the data empha- 
size that the observed magnitude-frequency relation- 
ships of landslides skew volumetric erosion towards 
the larger, more infrequent, landslides [14,15], This 
is reflected in all nine simulations, where the median 
50% of the CRN-derived and volumetric erosion rates 
generally fall below the total imposed erosion rate in 
the simulation. The outer statistical bounds of the 
CRN-derived and volumetric erosion rates are sub- 
stantially higher than the imposed erosion rates, re- 
flecting the the infrequency, but importance, of these 
large events in controlling erosion rates in landslide 
dominated catchments. 

The damping and averaging effect of erosion rates 
derived from CRNs in sediment are also illustrated 
in the statistical spread of volumetric- versus CRN- 
derived erosion rates at all catchment scales, sediment 
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Figure 5: Box-and-whisker plots of CRN-derived erosion rates and volumetric erosion rates for nine separate simulations with three combinations of sediment detachment and total 
erosion rates. Upper solid line on each graph indicates total erosion rate ; lower solid line indicates sediment detachment rate. See text for detailed discussion. 
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detachment rates, and landslide rates, where the dis- 
tribution of CRN-derived erosion rates is tighter than 
contemporaneous volumetric erosion rates. This ef- 
fect is particularly notable at small to moderate catch- 
ment scales in regions with low rates of erosion by 
landsliding. In particular, at 1 mm/yr total erosion 
rate, the CRN-derived rates have a 50% smaller dis- 
tribution than the volumetric erosion rates. 

The effect of spatial averaging on the spread of 
volumetric- and CRN-derived erosion rates is also 
highlighted in our results (Fig. 5). The statistical 
spread of both volumetric erosion rates increases with 
increasing catchment size, presumably related to the 
increased likelihood of experiencing a large mass- 
wasting at greater catchment area, until a threshold is 
reached, at which point the catchment becomes large 
enough to adequately average large landsliding events 
and areas unaffected by mass wasting, and the spread 
of the data drops significantly. This spatial scale in 
our simulations appears to occur between fourth- and 
fifth-order catchments (a jump from ~ 8 km 2 to ~ 70 
km 2 ). The exception to this rule is at low total erosion 
rates, where landslides appear to occur infrequently 
enough at all catchment scales to never significantly 
increase the statistical distribution of erosion rates at 
third- through fifth-order catchments. 

Finally, given the damping and spatial averaging 
affects described above, it is heartening to note that 
the median 50% of observations are more or less 
consistent between volumetric and CRN-derived ero- 
sion rates (Fig. 5). At small catchment scales, 
CRN-derived rates are typically higher than volumet- 
ric rates (this is particularly clear at low total ero- 
sion rates), but converge at larger catchment scales 
(Fig. 6). The catchment size at which this con- 
vergence occurs varies, and decreases with increased 
rates of erosion. As a general rule of thumb, how- 
ever, it would appear that CRN-derived rates of ero- 
sion from sediments are statistically representative of 
volumetric rates in our numerical simulation at third- 
, or at most, fourth-order catchment scales. This 
observation indicates that CRN-derived erosion rates 
from sediments in landslide-dominated catchments 
may in fact be useful for looking at basin-wide ero- 
sion. Even modest-scale drainages, such as San An- 
tonio Creek, have several fourth-order catchments. 
Measuring CRN-derived erosion rates from each of 
these catchments should allow the identification of 





Volumetric Erosion Rate (mm/yr) 


Figure 6. Scatter plots of median cosmogenic versus median volu- 
metric erosion rates for (from top to bottom) total erosion rates of 1 . 
5, and 10 mm/yr). All plots indicate that cosmogenic erosion rates 
are higher than volumetric erosion rates, and illustrate the conver- 
gence of the two rates with increasing sediment detachment rate. 
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outlying rates that may be due to recent mass wasting 
events. Such outlier rates would approximate the re- 
cent volumetric erosion rates from the basin (Fig. 5). 
Nonetheless, even fourth-order catchments are likely 
to underpredict the overall erosion rate by 20 - 40%, 
due to the impact of rare, but unusually large, land- 
slide events. 

3.2. Simulated bedrock erosion rates 

In addition to tracking the CRN concentration at 
46 catchments within the simulation, we also calcu- 
lated the surface CRN concentration at 100 points 
randomly distributed throughout the model space, 
and analyzed such concentrations for erosion rates as 
though bedrock samples were collected from each of 
these points (eqn. 4). The statistical distributions 
of these erosion rates are derived from analysis of 
100 random points, sampled 100 times at 1,000 year 
intervals, over a 100.000 year model run, such that 
each box-and-whisker in the plot represents the sta- 
tistical distribution of some 10,000 independently de- 
termined erosion rates (Fig. 7). 

A first-order observation is that the spread of ero- 
sion rates increases with increasing rates of total ero- 
sion (and thus, with increasing rates of landsliding). 
This almost certainly reflects the increased likelihood 
of any of the randomly sampled points in the model 
to be affected by landsliding with the increasing fre- 
quency of events. A second important observation is 
that at low rates of total erosion (and landsliding), a! 
though sediment detachment encompasses, at most, 
30% of total erosion, the median CRN-derived ero- 
sion rate is almost identical to the sediment detach- 
ment rate. That is, for any given point in the land- 
scape. the likelihood of sampling a point that has been 
recently enough effected by mass wasting to alter the 
CRN concentration is virtually negligible. As the to- 
tal rate of erosion increases, the median erosion rates 
increases above the background sediment detachment 
rate, yet fall well below the total erosion rate over the 
landscape. In contrast to the distribution of sediment- 
derived CRN erosion rates (Fig. 5), the upper extent 
of the ‘whisker’ (3 inter-quartile ranges beyond the 
median) never exceeds the imposed total erosion rate. 
Although the CRN-derived erosion rates do not fol- 
low a Gaussian distribution, one may draw the anal- 
ogy that less than 99% of the ‘bedrock’ samples will 
fail to accurately reflect the total erosion rate over the 


landscape. (The highest actual percentage is 0.25% 
for a total erosion rate of 10 mm/yr, and a sediment 
detachment rate of 0.01 mm/yr, where 25 samples 
out of 10,000 fell within 2 mm/yr of the total erosion 
rate). 

These results suggest that sampling bedrock expo- 
sures in basins dominated by mass wasting may pro- 
vide an upper bound on sediment detachment rates 
across the basin, but are inadequate as a means to de- 
rive overall rates of erosion from the combined effects 
of sediment detachment and landsliding. 

33. Response of CRN-derived erosion rates to 
changes in rates of mass wasting processes 

In addition to tracking the CRN concentrations in 
sediment from catchments, and at individual points, 
the simulation also calculates the mean CRN concen- 
tration across the landscape at each point in the sim- 
ulation. Whereas the calculation of CRN and volu- 
metric erosion rates is only performed while the sim- 
ulation is in an erosional ‘steady -state’ for the im- 
posed landslide and sediment detachment rates, the 
mean CRN concentration is calculated throughout the 
model run to gather an estimate of the response time 
of CRNs in the landscape to changes in erosional 
boundary conditions (Fig. 8). 

These plots record the mean 10 Be concentration 
over the landscape, beginning at the initialization of 
the model, where the 10 Be concentration at each cell 
in the simulation is analytically solved for the im- 
posed sediment detachment rate and the scaled pro- 
duction factor at each cell. Subsequently, landslides 
are populated across the model landscape, and the 
mean concentration of 10 Be begins to decrease at a 
rate controlled by the sediment detachment rate and 
the total rate of landsliding imposed on the model. 
After a period when mean 10 Be concentration steady- 
state is achieved during sediment detachment and 
landsliding, the landslides are eliminated from the 
simulation, and the mean 10 Be concentration in the 
landscape increases (Fig. 8A). 

Tracking of specific points emphasizes the effects 
of changes in sediment detachment and landsliding 
rates on the response time of mean 10 Be concentra- 
tion (Fig. 8B). The response time of the landscape 
to achieve a new equilibrium mean 10 Be concentra- 
tion is a function of both the sediment detachment 
rate and the landsliding rate. For a given total ero- 
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Sediment detachment rate (mm/yr) 


Figure 7. Distribution of cosmogenically determined erosion rates (shaded boxes) from individual sample points for each model run. Total 
erosion rates for each set of runs are indicated by solid lines; sediment detachment rates are drawn across the lower portion of the graph 
in dashed lines. As landslide erosion rates increase, the spread and median of cosmogenically derived erosion rates also increase. Median 
cosmogenically determined erosion rates at these points, however, remain well below total erosion rates for each model run. 


sion or sediment detachment rate, the time for 10 Be 
to reach a new equilibrium decreases with an increase 
in the other rate. Increasing the sediment detachment 
rate by a factor of ten decreases the time necessary 
for 10 Be to reach a new equilibrium concentration 
(squares on Fig. 8B); likewise, increasing the rate of 
landsliding by a factor of 5 causes also yields a de- 
crease in the response time (circles on Fig. 8B). On 
the other hand, the response time for re-equilibration 
of mean 10 Be concentration at the cessation of land- 
sliding is a function solely of sediment detachment 
rate. The time necessary for a landscape to recover 
95% of its pre-landslide mean 10 Be concentration is 
as little as 2,000 yrs, at sediment detachment rates of 
0.3 mm/yr, and as great as 55,000 to 75,000 years at 
sediment detachment rates of 0.01 mm/yr (triangles 
on Fig. 8B). 

These results indicate that in rapidly eroding land- 
scapes, changes in landsliding rates over relatively 
short time periods (say a few thousand years or less) 
may generate transients in CRN concentrations last- 
ing tens of thousands of years. If landsliding is a 
more prevalent mechanism of erosion under certain 
climatic regimes (e.g. during interglacial periods; 
[21,?]), then the response time of mean 10 Be concen- 
tration to changes in landsliding rate must be consid- 


ered when the frequency of these rate changes is of 
order the 10 Be response time of tens of thousands of 
years. 

3.4. Comparison of simulation results of CRN ero- 
sion rates to CRN erosion rates from the 
Nepalese Himalaya 

The Khudi River in Nepal, a tributary to the 
Marsyandi, has been the focus of a multi-disciplinary 
study of geomorphic and geodynamic coupling in the 
Flimalaya. As part of this study, erosion rates over 
this catchment have been assessed using a variety of 
techniques, including long-term erosion rates from 
low-temperature thermochronometers (> 2-5 mm/yr; 
1 22]), and present-day erosion rates from stream- 
sediment suspended-load determinations (~ 3 mm/yr; 
E. Gabet, pers. comm.). In addition, 10 Be CRN ero- 
sion rate determinations were made throughout the 
catchment, both from bedrock exposures and 0 th - 
order (< 0.01 km 2 ) catchments, and from the mouth 
of the Khudi River where it joins the Marsyandi 
(Heimsath, in prep.). The measured CRN erosion 
rates were calculated using the same production rates 
and scaling factors as in the model, and have an av- 
erage error of 10%. Here we compare the statistical 
distribution of 56 10 Be CRN erosion rate determina- 
tions from bedrock exposures and 0 th -order catch- 
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Figure 8. Mean concentration of I0 Be across the simulation, showing the effect of initiation and cessation of landslides on 1 °Be concentration 
at several different rates of sediment detachment and landsliding. A) Mean I0 Be concentration through time for all model runs; initiation and 
cessation of landsliding indicated. B) Symbols highlighting variation in response time to changing boundary conditions, see text for discussion. 
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merits with predicted statistical distributions of CRN 
erosion rates determined from our simulation. Ad- 
ditionally, we calculate distributions of CRN erosion 
rates predicted at the mouth of the Khudi River, and 
compare these with the basin-wide CRN erosion rates 
determined for the Khudi catchment. 

The simulation of landsliding, cosmogenic produc- 
tion, and erosion was performed as described for the 
theoretical study in the San Antonio Creek watershed, 
with only a few modifications. First, the highest reso- 
lution DEM available for Nepal has roughly 90-m cell 
spacing, as opposed to the 30-m spacing available for 
San Antonio Creek. The model cell size was adjusted 
accordingly, and the smallest landslide populated in 
the model was also adjusted to have a radius of 90 
meters. Second, available erosion rate data from ther- 
mochronologic and sediment-load studies were used 
to establish a range of landsliding rates to input to the 
model of the Khudi catchment. A sediment detach- 
ment rate of 0.15 mm/yr was assumed for the initial 
model runs, and landsliding rates of 2.85, 3.35, and 
3.85 mm/yr (for total erosion rates of 3, 3.5, and 4 
mm/yr) were selected. The best-fit run was then re- 
analyzed with varying sediment detachment rates to 
derive a statistical distribution of bedrock and small- 
order catchment CRN erosion rates for comparison 
with existing CRN data. Finally, the size of the Oth- 
order catchments is of order the size of the model 
cell spacing, so these small catchments were treated 
as points within this model. 

3.4.1. Basin-wide CRN erosion rate in the Khndi 
catchment 

The distribution of basin-wide CRN-erosion rates 
derived from the model were compared with a 10 Be- 
derived erosion rate measured on sediment deposited 
at the mouth of the Khudi catchment. The measured 
CRN erosion rate of 3.4 mm/yr matches the median 
of the CRN erosion rates for the model run with total 
erosion rate of 3.5 mm/yr (0.15 mm/yr sediment de- 
tachment and 3.35 mm/yr landslides; Fig. 9A). This 
result supports the theoretical determination that at 
high landslide erosion rates, the effects of landslid- 
ing on CRN erosion rates are spatially averaged over 
large catchments, and that CRN methods are an effec- 
tive means of assessing average erosion rates. Further, 
the median values of each of these three model runs 
does not overlap the 1 st through 3 rd quartiles of any 


other run (Fig. 9A), indicating that the medians of 
these model runs statistically differ [23], The results 
of this simulation, in concert with the actual CRN ero- 
sion rate, allow a determination of the erosion rate in 
the Khudi catchment of 3 ± 0.5 mm/yr, tightening the 
constraints on the erosion rate as derived from ther- 
mochronologic and sediment-load data. 

3.4.2. Bedrock and small-order catchment CRN 
erosion rates in the Khudi catchment 

Using the basin-wide total erosion rate of 3.5 
mm/yr, four additional model runs were executed at 
increasing rates of sediment detachment (0.01, 0.05, 
0.10, and 0.15 mm/yr) to compare the simulated dis- 
tributions of CRN erosion rates with the distribution 
of CRN erosion rates determined from 56 10 Be ero- 
sion rates from bedrock samples and 0 th -order catch- 
ments (Heimsath, in prep.; Fig. 9B). In the model 
space, the same 56 locations that were actually sam- 
pled were used to determine a theoretical distribution 
of bedrock erosion rates. The results of this compar- 
ison indicate that, as hypothesized from the theoreti- 
cal results, sampling bedrock and small area drainage 
basins for CRN concentrations in regions with sig- 
nificant erosion by bedrock landsliding is not an ef- 
fective approach to determining basin-wide erosion 
rates. The data set analyzed in the Khudi drainage is 
one of the larger CRN erosion rate studies undertaken, 
yet the number of samples collected, and the percent- 
age of the landscape that they represent, is inadequate 
to provide a meaningful representation of basin-wide 
erosional processes. Although thermochronometric, 
cosmogenic, and sedimentologic techniques all indi- 
cate erosion rates in the Khudi of 2-5 mm/yr, the me- 
dian CRN erosion rate from the Khudi bedrock sam- 
ples is ~ 0.1 mm/yr. while the median of the simula- 
tion runs is ~ 0.25 mm/yr (Fig. 9B). It is also interest- 
ing to note that the median erosion rates from the four 
model runs are very similar, despite an order of mag- 
nitude difference in sediment detachment rate among 
them, and they differ significantly from the median 
erosion rate determined on the actual bedrock sam- 
ples. This difference can likely be ascribed to two 
sources, both due to the fact that the actual bedrock 
samples from the Khudi catchment almost all came 
from ridge crests, or small draws that head at ridge 
crests. The first source of the observed difference 
is that landslides in the model space are placed ran- 
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Figure 9. A) Distribution of model basin-wide CRN-derived erosion rates from the mouth of the Khudi drainage for erosion rates of 3.0, 3.5, 
and 4.0 mm/yr. The thick gray line across the center of the figure is the measured 10 Be derived erosion rate of 3.4 mm/yr from the mouth 
of the Khudi River (A. Heimsath, unpub. data). B) Comparison of simulated distributions of CRN-derived erosion rates from 56 sample 
locations (white boxes) where actual CRN-erosion rates were measured (gray box: A. Heimsath. unpub. data) for bedrock samples in O th -order 
catchments in the Khudi drainage. Median CRN erosion rates for both the simulated data and the actual CRN data are an order of magnitude 
less than the basin-wide erosion rate (A). 


domly, such that a ridge crest is as likely to be af- 
fected by landslides as a hillslope. In reality, this is 
likely to not be the case. Landslide mapping indicates 
that although landslides are preferentially located in 
areas of steep slopes near the headwaters of drainage 
systems, such slides rarely breach the drainage di- 
vide and lower the interfluve (e.g. [24]). Lowering 
of drainage divides most plausibly occurs during in- 
frequent events when the topographic slope below' the 
ridge crest has been over steepened by repeated land- 
sliding, and is thus this process is likely to be highly 
undersampled. The second effect is sampling bias. 
The CRN samples were specifically collected from 
bedrock outcrops along the crest that were judged to 
not have been recently affected by landsliding. As 
such, the spread of CRN ages would be expected to 
be smaller than a truly random set of ridge crest sam- 
ples, and the median erosion rate would, as a result, 
be lower. This difference notwithstanding, the results 
of this comparison still serve to emphasize the impor- 
tance of considering the effects of spatial and tem- 
poral averaging in collecting material for CRN deter- 
minations of erosion rates in actively landsliding re- 


gions. 

4. Discussion and Conclusions 

We have presented a numerical simulation for mod 
eling the production, decay, and distribution of cos- 
mogenic nuclides on a landscape, and their removal 
through sediment detachment and mass wasting pro- 
cesses. Although this simulation was developed to 
model the effects of landsliding on CRN-derived ero- 
sion rates, the framework of the simulation could 
be adapted to model any number of factors that ef- 
fect erosion rates derived from CRNs in sediments, 
including spatial variations in lithology and mineral 
content, ice cover, annual snowfall, and recent effects 
of glaciation. In the past, careful researchers trying 
to exploit CRNs to obtain erosion rates have com- 
monly restricted their sampling to small, unglaciated 
catchments with uniform lithologies and slow' ero- 
sion rates. Variations in CRN production due to to- 
pographic shielding, slope, lithology, snow and ice 
cover, and glacial history have generally been ig- 
nored. With the ability to numerically predict the 
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effects of such variations, studies can be expanded 
into larger catchments with higher erosion rates, di- 
verse sediment-production processes, and spatial het- 
erogeneities in CRN concentrations and events that 
can re-set the cosmogenic clock. 

Our landsliding model is underpinned by the ob- 
served frequency-magnitude relationships of land- 
slides in two mountain ranges. These happen to 
yield nearly identical exponential scaling factors (0 = 
1.1) that imply that large, infrequent landslides dom- 
inate the total sediment flux. The observation that 
implementation of our frequency-magnitude based 
landsliding module yields volumetric rates that os- 
cillate around the expected value suggests that the 
model succeeds in mimicking a natural process. This 
frequency-magnitude relationship and the appropriate 
scaling factor need to be verified in other mountain 
ranges before they are routinely applied. In addition, 
we make the implicit assumption that all of the sedi- 
ment generated by a landslide is delivered and homo- 
geneously mixed with sediment detachment-derived 
sediment within a single model time step (100 years). 
We make no attempt to model sediment storage on 
hillslopes or within fluvial systems, or to model parti- 
cle size fractionation between sediment detachment- 
derived and landslide-derived debris. 

Our specific evaluation of the effects of bedrock 
landsliding on erosion rates derived from CRNs us- 
ing this simulation yielded several results that have 
been previously described, particularly that stochas- 
tic processes, dominated by large rare events, are 
difficult to measure using basin-averaging sediment- 
sampling techniques, because the large, rare events, 
will often not be represented in the sample popula- 
tion (e.g. [?]). However, our simulations also indi- 
cate that the median CRN-determined erosion rates 
are representative of the volumetric erosion rates de- 
rived from the same catchments. This suggests that 
sampling multiple, similarly sized catchments, even 
in active, landslide-dominated mountain belts, offers 
a significant likelihood of yielding several samples 
with consistent CRN erosion rates. These rates are 
likely to be representative of the recent erosion within 
those catchments, although it must be recognized that 
such results will typically be lower than the long- 
term average that samples the large, rare events. As 
long as spatial variations in production can be ade- 
quately accounted for, as is done within our model- 


ing environment, then larger catchments will always 
yield a better approximation of long-term erosion rate 
in landslide-dominated terrains than smaller catch- 
ments. Although a specific relationship between the 
catchment size necessary to spatially average CRN 
samples and erosion rate is difficult to derive, based 
on our modeling results a general rule of thumb ap- 
pears to be 

Aw 9 = ^ (7) 

where A avg is the area needed to average the vari- 
ability in CRN concentration, and E is the estimated 
erosion rate over the catchment. 

Sampling from bedrock outcrops to measure aver- 
age erosion rates in landslide-dominated catchments, 
however, is unlikely to be a useful exercise. At the low 
end of landslide erosion rates, such samples will faith- 
fully yield the sediment detachment rate on the land- 
scape, but at increased rates of landslide erosion, such 
samples will only yield a very rough upper bound on 
the sediment detachment rate, and, more likely, will 
be uninterpretable in the context of the spatially aver- 
aged erosion rate over the study area. 

Finally, the response time of 10 Be concentrations 
over the landscape to changes in rates of erosional 
processes is thousands to tens of thousands of years. 
Regions that have undergone recent changes in rates 
of erosion may yield CRN-derived erosion rates that 
reflect some intermediate rate between the previous 
and current erosion rates during the re-establishment 
of ln Be equilibrium over the landscape. 
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Appendix A 
Governing Equations 

This Appendix details the derivations of three spe- 
cific portions of the landslide model, calculating the 
rate of landsliding, k, given a mean landslide erosion 
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rate, Ei s \ population of landslides in model space, and 
determining mean cosmogenic concentration per vol- 
ume of material removed from the model. 

Rate of Landsliding 

We assume that landslides follow a power law fre- 
quency magnitude distribution [7,14]. For such a dis- 
tribution, the cumulative frequency of landslides can 
be written as 

tia>a. = * (A a /A T )~ e A r , (A8) 

where tia>a. is the number of landslides greater than 
area A s that occur in a given year, A r is a specified 
reference area, k is the rate of landsliding, and 8 is 
the power-law exponent of the frequency-magnitude 
distribution [7, see eqn. 1], If A r is taken to be 1 km 2 , 
then from Eqn. A8, the total number of landslides per 
year, n, over the reference area can be written as 

n T = f nA~ p dA, (A9) 

“ Amin 

= K | A-efr", (A 10) 

= *{KL-Ki t)- (AH) 

where A min is the area of the smallest landslide to oc- 
cur in A r (in the case of the simulation, A min is equal 
to the model cell size) and A max is the largest slide 
area, constrained by local topographic relief. Since 
A m i n and A max can be constrained either empiri- 
cally or based on physical characteristics of the modei 
space, and /? has been shown to be ~ 1 over a large 
range of landslide erosion rates [7.14,15], the volume 
of material removed by landsliding must be controlled 
by the rate of landsliding, k. Because we would prefer 
to prescribe a rate of erosion due to landslides, Ei s for 
the model, we need to solve for k in terms of Ei s . To 
do this, we start by determining the number of slides, 
tia „, of a given area, A s per year over the reference 
area, A r : 

tia , = (A 12) 

Given this, the volume of erosion in any year due to 
landslides of area A s is equal to the number of land- 
slides of area A s multiplied by the volume of land- 
slides of area A„, 

V A , = tia. • V s (A 13) 

= nPA-'-e-V.. (A 14) 


In this simulation, we have opted to model the land- 
slides with parabolic cross sections and a linear re- 
lationship between maximum slide depth and width 
[25,7]. The scaling between landslide area and depth 
is defined by a scaling factor, £, with an empirically 
determined value of ~0.05 [7]. Thus, the depth of 
the landslide, d, at any radial distance from the land- 
slide center is a function of the maximum slide depth, 
the radial distance from the landslide center, r, and a 
constant, C, 

d = d max — Cr 2 . (A15) 

We can solve for C at the outer edge of the landslide, 
r m «i, where d = 0, 


For a slide of a given area A s , then, the maximum 
landslide depth, d max and radius, r max , are given by 


dmax — £ x/ As 


Substitute these values into Eqn. A16 to solve for C 
in terms of area, A s : 


and substitute Eqns. A17 and A19 into Eqn. A15 to 
solve for d as a function of A s , 

d = £-/a 7 7==r 2 . (A20) 

VAs 

The volume, V s , of a landslide of area A s can then be 
calculated by integrating over cylindrical shells from 
V — 0 tO T ma x , 


r / SyJ~A s - 
J 0 


=r 2 • rdr (A21) 


0 | 2 ^4 

= 2 tt £ r 7=r 

I 2 4VC47 0 

9 (aT ^ /2 \ 

= 2n£ {-2^--^r) 

= Ui'* 

2 8 
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Eqn. A24 can be substituted back into Eqn. A 14 to 
yield the volume of material, Va., removed from the 
reference area, A T per year by slides of area A s \ 

V A . = \nep AW*-®. (A25) 

The total volume of erosion per year from reference 
area A r , then, is the sum of the erosion due to land- 
slides of all sizes that occur in A r , from A m i n to 
Amax > 




1 r^max 

Ike/ 3 [ A\] /2 - B) dA (A26) 

2 JA mi „ 

«A27, 


The total erosion rate per year due to landslides, then, 
can be determined by dividing the volume of material 
eroded per year by landslides over A r by the area of 
A r , 



(A28) 


and, since we have defined A r to be a unit area (1 
km 2 ). 


Eu = V r (A29) 

Thus, Eqn. A27 can be solved for k in terms of £) s to 
yield 


* = -/ 'I, , )V (A30) 

Substituting Eqn. A30 into Eqn. All yields the to- 
tal number of landslides, n r , per year over the refer- 
ence area, A r in terms of known or prescribed values 

Amin, Am ax , &■> Und El s ', 


Tl r 


E u (3-2/3) (A-^-A-f n ) 


(A31) 


The total number of landslides, n; s , then, in the simu- 
lation for a given time step, t, over the entire simula- 
tion is 


n ls = t x n T x 4=2. (A32) 

A r 

where A S j m is the area, in km 2 , of the simulation. 


Landslide Population and Distribution 

With the number of landslides per time step deter- 
mined, the simulation is populated. The position each 
landslide is specified by a randomly generated x,y co- 
ordinate pair. The size of a given landslide is derived 
from the landslide frequency-magnitude relationship 
(Eqn. A8). The probability of a landslide with area 
A„ occurring is 

Pa. = « A~ 0 (A33) 

therefore, randomly generated numbers mapped lin- 
early onto the range Pa„ 4 „ to P Ama , can be used 
to create a population of nu landslides that fit the 
frequency-magnitude distribution nA~ 0 . 

Cosmogenic Nuclide Concentration in Eroded Ma- 
terial 

For each model step, the total depth of material re- 
moved from a given cell, D, is the sum of material 
eroded by sediment detachment, E s and the material 
eroded by landslides, Ei s . Using the surface con- 
centration of a cosmogenic nuclide at the end of the 
model step, A r ^°Be, the average nuclide concentration 
in the volume of eroded material, AT^Be can be cal- 
culated by integrating over the total depth of eroded 
material the concentration of the nuclide as a function 
of e-folding depth, A, 


AT)°Be • p , 

[ e~ z ^ A) dz 

(A34) 

O 

-D 


Nj°Be ■ p ■ 

A L e -z(p/A)|° 

(A35) 


p \ 1 -D 


NfBe • A ( 

r e D(p/A) _ 

(A36) 
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